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We formalize the Hubble slow-roll formalism for inflationary dynamics in Randall-Sundrum Type 
II braneworld cosmologies, defining Hubble slow-roll parameters which can be used along with the 
Hamilton-Jacobi formalism. Focussing on the high-energy limit, we use these to calculate the exact 
power spectrum for power-law inflation, and then perturb around this solution to derive the higher- 
order expression for the density perturbations (sometimes called the Stewart-Lyth correction) of 
slow-roll braneworld models. Finally we apply our result to specific examples of potentials to 
calculate the correction to the amplitude of the power spectrum, and compare it with the standard 
cosmology. We find that the amplitude is not changed significantly by the higher-order correction. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



New ideas in fundamental physics have suggested the 
possibility that our observable Universe lies on a 'brane' 
within a higher-dimensional bulk space-time, and this 
idea may have serious ramifications for early Universe 
cosmology pj. The most studied scenario, viewed as a 
toy model for more elaborate proposals, is the Randall- 
Sundrum Type II (RS-II) model where there is a sin- 
gle brane living in an anti de Sitter bulk 0. In that 
case, the Friedmann equation is modified at high energy, 
potentially with significant consequences for our under- 
standing of inflation and the perturbations it generates. 

The most powerful tool for studying inflationary dy- 
namics is the slow-roll approximation which is ex- 
pected to be extremely accurate for models capable of 
matching current observations such as those of WMAP 
4]. This enables considerable analytical progress to be 
made, and predictions for the perturbations generated 
by inflation are most conveniently expressed in terms 
of slow-roll parameters which measure the accuracy of 
the approximation Such an approach was used by 
Maartens et al. [|| to describe the outcome of inflation in 
the RS-II braneworld scenario. 

In that paper, the version of the slow-roll approxima- 
tion used was based on derivatives of the potential driv- 
ing inflation. In this paper we formulate the approxi- 
mation in terms of derivatives of the Hubble parameter, 
which in the standard cosmology is an important tool 
for developing high-accuracy predictions for density per- 
turbations. We take advantage of our formalism to carry 
out similar calculations for the high-energy limit of RS-II 
inflation models, deriving an exact result for power-law 
inflation and calculating the higher-order slow-roll cor- 
rection to the perturbation amplitude (often called the 
Stewart-Lyth correction, as they carried out the original 
calculation for the standard cosmology 0). We find that 
the correction is of a similar size to that in the standard 
cosmology (i.e. negligible in many inflationary models), 
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and calculate it for some sample inflation models. 



II. SLOW-ROLL FORMALISM 

We follow the notation set down by Liddle and Taylor 
@. In the Randall-Sundrum Type II model @ the Fried- 
mann equation receives an additional term quadratic in 
the density 0- The Hubble parameter H is related to 
the energy density p by 



H 2 = 



8tt 
3M| 



P 



(1) 



where M4 is the four-dimensional Planck mass and A is 
the brane tension. We have set the four-dimensional cos- 
mological constant to zero, and assumed that inflation 
rapidly makes any dark radiation term negligible. This 
reduces to the usual Friedmann equation for p <C A. We 
assume that the scalar field is confined to the brane, so 
it obeys the usual equation 



3H<j) = -V' . 



(2) 



where prime indicates derivative with respect to <j), and 
dot a derivative with respect to time. Its energy density 
is p = V + 2 /2. 

For the standard cosmology, the Hubble slow-roll for- 
malism was set down in detail in Ref. 10] . The first two 
parameters are defined as 



£h = 3 



V + c/) 2 /2 



3H<j) 



Mj H' 2 
17 IP 



17 IT 



(3) 



(4) 



If these parameters are much less than one, they allow the 
neglect of the <f> term in the Friedmann equation, and the 
ip term in the scalar wave equation In addition, the 



2 



condition for inflation, a > is conveniently expressed as 
e H < I- 

We seek a generalization of those parameters suitable 
for use with the braneworld equations. We note that 
there is some arbitrariness to the definitions, at least in 
terms of the constant prefactor. However this can be re- 
moved by imposing the requirements that the parameters 
en and ryn continue to correspond to the conditions en- 
abling neglect of terms in the Friedmann and fluid equa- 
tions respectively, that eu < 1 remains the condition for 
inflation, and that in the slow-roll limit the density per- 
turbation spectral index takes the same form as in the 
standard cosmology. 

We will make use of an approach developed by Hawkins 
and Lidsey (Hill, 

who devised a formalism for brane- 
world inflation with many of the properties of the 
Hamilton Jacobi approach used in the standard cosmol- 
ogy |l3j| . They define a quantity y (</>), which is to play a 
similar role to H ((f)) in the standard cosmology, by 



P/2A 
1 + p/2\ 



The inverse relation is 



P ; 



2Xy 2 



1-y* 



(5) 
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0, while at high energies 



At low energies p <C A, y 2 
p> A, y 2 f. 

The Friedmann equation can be written as 



H(y) 



16?^ 1/2 .// 
3Mf" 



1 



y 



(7) 



From the scalar wave equation and Eq. © one obtains 

y' 



(xm^ 1/2 



V 3tt 



i-y 2 



(8) 



Taking the derivative with respect to the field in the 
Friedmann equation, one finds 



H' = 



/16vrA\ 1/2 


i + y 2 


{SMI) 


L(l-y 2 ) 2 J 



and using Eq. 



H' 



4tt (l+y 2 ) 



(9) 



(10) 



M 2 (l~y*r- 
In Ref. 12] , Hawkins and Lidsey define two parameters 



Mjy^ 
Ait y 2 
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(11) 



by analogy to the standard cosmology Hubble slow-roll 
parameters. These parameters prove useful in analyzing 
the exact dynamics of braneworld inflationary models. 
However, their smallness (as compared to unity) does 



not precisely correspond to the ability to neglect terms 
in the Friedmann and scalar wave equations, and this 
means they are not ideal for the purpose of analyzing 
perturbation generation. We therefore define new Hubble 
slow-roll parameters for the braneworld, which do have 
this property. 

To do this, we define the parameters as ratios of terms 
in the Friedmann and wave equations, following Eqs. © 
and igj: 



£H = C(y) 



V + <t> 2 /2 



; m = -D(y) 



(12) 



At this stage we have allowed the 'constant' prefactors 
C{y) and D(y) to depend on y; we will next show how to 
fix them using the requirements that en < 1 corresponds 
to inflation, and that in the slow-roll limit the density 
perturbation spectral index takes its usual form. In or- 
der for these parameters to correspond to the ability to 
neglect terms, those prefactors should always be of order 
unity, and we will soon see that they are. 

From the definition of en, using Eqs. and fTTijl . one 
can find 



C(y) (XM 2 
3 V 3?r 



1/2 



y 



H'- 



{l + y 2 ) 2 H 3 



(13) 



The coefficient C(y) can be determined by demanding 
that ci = ■<=>■ en = 1; taking the derivative of the 
Friedmann equation with respect to time gives 



a 
a 



16nXy 2 



Mf (1 - y 2 f 



-q + y 2 ) 
C(y) 



(14) 



and so we require C(y) = 3(1 + y 2 )- Its value ranges from 
3 in the low-energy limit to 6 in the high-energy limit. 
Our definition for en is therefore 



1/2 



H'- 



( XM 2 
\ 3vr J (l + y 2 ) H 3 



(15) 



In the low-energy limit, Eq. I|15f) becomes the usual ex- 
pression Eq.J!]J). 

For t?h we apply Eqs. (Q and l|l(J|l to obtain 



D{y) (XM 2 
3 ^ 3tt 



1/2 



H" 



4y3 



H 
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(l + y 2 )H 2 (l + y 2 ) 3 H 3 



There are several ways one could aim to fix the constant 
D(y) to establish a unique definition of rju- We choose 
to do so such that the slow-roll expression for the density 
perturbation spectral index takes its usual form, namely 
n = 1 — 4en + 2t;h, as was done for the potential slow-roll 
parameters in Ref. The slow-roll expression for the 
perturbation amplitude is 



-,1/2 



H 2 
2tt</> 



(17) 



iH 
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which in terms of those variables is exactly the usual 
result. The spectral index is defined as 



n-l = 



dlnV-R 



din fc 

and using Eqs. (JJJ and (|10fl one obtains 



(18) 



n-1 = -4 



AMf 
3tt 



1/2 



y 



H' 



( AM, 



V 3tt 



1/2 



(1 + y 2 ) H 3 

y H" 



(19) 



4y 3 



(l + y 2 ) H 2 (1 + y 2 ) 3 i? 3 



This agrees with the usual expression n = 1 — 4en + 2?7h 
provided £>(?/) is set to 3 for all regimes. That gives us 
our definition 



XMj 
3tt 



1/2 r 



IS 



H" 



4y 3 



(l + y 2 ) H 2 {l + y 2 ) : 



H<21 
IP 



(20) 



This too reduces to the usual expression, Eq. (@J, in the 
low-energy regime. 

The two definitions Eqs. (p~5f) and (ffilj) define Hubble 
slow- roll parameters valid in all regimes of RS-II brane in- 
flation, generalizing the usual ones while preserving many 
key results: they give the conditions for neglecting terms 
in the Fricdmann and fluid equations, en < 1 corresponds 
to an inflationary expansion, and the slow-roll spectral 
index formula is always n = 1 — 4ch + 2?7h- 



III. EXACT AND HIGHER-ORDER 
PERTURBATIONS IN THE HIGH-ENERGY 
REGIME 

In this section we exploit the formalism of the previ- 
ous section to make accurate calculations of the density 
perturbations. Throughout this section we will restrict 
ourselves to the high-energy regime, obtained by taking 
y — > 1, where our slow-roll parameters can be written 



Mf H' 2 
4tT IP 



17 1p 



(21) 



where Ms = (AirX/S^^M^ 3 is the five-dimensional 
Planck mass. The high-energy versions of Eqs. (|10[) and 
are 



and 



H 



Ml H' 
17 ' 



47T 

3MI 1 



(22) 



(23) 



All these expressions could have been obtained directly 
for the high-energy regime using the same criteria we set 
down for the general case in the previous section. 



We will consider one case where the perturbations can 
be obtained exactly (namely power-law inflation, though 
in this case the corresponding potential is not exponen- 
tial), and then carry out our main calculation which is to 
compute the correction to the density perturbation am- 
plitude from next-order in slow-roll. This type of cal- 
culation was first performed by Stewart and Lyth for 
the standard cosmology |7J, and is often known as the 
Stewart-Lyth correction. We will compute its equivalent 
for the density perturbations in the high-energy regime 
of the RS-II model. 

We will calculate the perturbations using a formalism 
due to Mukhanov 01 ■ He defined a new variable u = aStfi 
(where Sep is defined in the spatially-flat gauge; a gauge- 
invariant definition can be made which includes a con- 
tribution from the curvature perturbation), and demon- 
strated that in linear perturbation theory its Fourier 
modes obey the wave equation 



d 2 u k 
dr 2 



k 2 --—} u k = Q. 

z dr z 



(24) 



Here r is the conformal time, and z = a<j)/H encodes 
all the relevant information about how the background 
is evolving. This equation has been fully derived only in 
the standard cosmology, and we stress that it is currently 
a conjecture that it can still be used in the braneworld 
context. This is because as yet no-one has been able to 
calculate the backreaction from five-dimensional gravity 
in order to assess whether it is different from the four- 
dimensional backreaction that the Mukhanov equation 
incorporates, other than in the large-scale limit where 
energy conservation alone is sufficient to ensure the per- 
turbations remain constant |lfij |. There is however sup- 
porting evidence for its use, because its derivation does 
not need the Friedmann equation. This is clear since the 
Mukhanov equation does not feature the gravitational 
constant G (though some derivations of it do appear to 
use the Friedmann equation in intermediate steps), and 
we have verified explicitly that the same equation does re- 
sult using the modified Friedmann equation in the deriva- 
tion as given in Ref. In other words, the Mukhanov 
equation tells us how scalar field perturbations develop in 
a given expanding background, without needing to know 
the physics responsible for determining the background 
evolution. This however falls short of being a full five- 
dimensional calculation as would be required to fully ver- 
ify the use of the equation, but as yet it is not known how 
to implement such a calculation. 

Having adopted the Mukhanov equation, the first step 
is to write it in terms of the slow-roll parameters, with a 
lengthy calculation yielding 



1 d 2 z 
z dr 2 



= 2a 2 H 2 



3 7 
1 + e H - -m - 2 £ h??h+ 



1 



Ml H'H 1 ' 



32tt 2 H 4 



(25) 



which is an exact relation. 
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A. Exact mode equation solution 

It is well known that the mode function can be solved 
exactly if the square-bracketted term in Eq. I|25|) is con- 
stant, which in the standard cosmology corresponds to 
power-law inflation from an exponential potential |17| . 
the calculation having first been performed by Lyth and 
Stewart • To discover if there is an analogous result 
for the RS-II scenario, we set en = l/p where p > 1 is 
a constant, and examine whether this makes the square 
bracket constant. 

Taking the derivative of en with respect to the field, 
one can write 



2ry H - eh] 



(26) 



As we have demanded en is constant, this implies 
?;h = l/2p. Similarly, differentiating rj^ gives 



Mf H'H' 
32ir 2 H 4 



{m + eh) 



H 
2H 



7 (»7h + 4) 



implying 



MgH'H'" 
32tt 2 # 4 



3 

v 



(27) 



(28) 



The square-bracket of Eq. Q25JI therefore indeed is con- 
stant, so the equation can be solved exactly. 

Before going on to do that, it is interesting to ask what 
potential gives this solution. Solving for the Hubble pa- 
rameter from the definition of en, and then substitut- 
ing into the Hamilton-Jacobi equation in the high-energy 
limit, namely 



M$H 



An 



3M| 



we find that the corresponding potential is 

lM 5 6 (6p-l) 



V{<f>) 



8 7T 2 2 



(29) 



(30) 



Instead of the exponential potential found in the stan- 
dard cosmology, we have an inverse power-law potential. 
Nevertheless, the corresponding expansion law a cx t p is 
power-law inflation as usual. 

Following Refs. 0,0,0], the conformal time for con- 
stant eh is given by 



1 



1 



aH 1 — en 



(31) 



and Eq. I|25[) with the values of the Hubble slow-roll pa- 
rameters gives 



1 d 2 z _ 1 8p 2 +2p- 1 
z dr 2 4 T 2 (p - l) 2 



(32) 



This allows one to write the Mukhanov equation Eq. 124|) 
as a Bessel-like one 



d 2 +fc2 ^ 1/4) 



Uk = 0, 



(33) 



with v — 3p/2(p— 1). The solution with the appropriate 
behaviour at small scales can be written as 

Uk {r) = | e *+i/^/ 2( _ r) i%(i) ( ^ T)] (34) 

where is the Hankel function of the first kind of 

order v. The asymptotic form of this equation once the 
mode is outside of the horizon is obtained by taking the 
limit k/aH ^ 



s i( J /-l/2)7r/2 2 i'-3/ 



2 m i 

T(3/2) V2fc 



(-fcr) 



-I/+1/2 



> (35) 



from which the corresponding form of the power spec- 
trum using |l9j 



yields 



V]l 2 {k) 



rp) W\ 



Ml 



(36) 



(37) 



k=aH 



B. Higher-order perturbation calculation 

We now perturb around the exact solution given above 
for small eh and rjn, following Refs. and [lSj . The 
expansion to lowest order of the conformal time gives 



(38) 



Applying this in Eq. (|25|l and truncating the expansion 
to first order, one arrives at another Bessel equation, now 
with v given by 



v ~ - + 2eh - »7h • 



(39) 



Note that the final three terms of Eq. 1251) do not affect 
the form of this expression. Then the solution Eq. I|35() 
can be used with the new form for v and the conformal 
time, expanding also the Gamma function and the other 
expressions to first-order, to obtain the final result 



,1/2 

n 



[1 - (2C + l)e H + C m ] 



2H 3 



Ml\H'\ 



(40) 



iH 



where C = -2 + In 2 + b ~ -0.73, with b the Eulcr- 
Mascheroni constant. The leading-order term, obtained 
by setting the square bracket to one, agrees with the 
result of Maartens et al. |(| in the high-energy limit. 
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While formally the correction term in the square- 
bracket looks exactly the same as in the standard cos- 
mology 0, E] , we should recall that the slow-roll param- 
eters which appear in it are generalizations of those in the 
standard cosmology, and could have been defined in dif- 
ferent ways. To get a feel for what this correction means, 
we need to evaluate it for some characteristic potentials, 
which we do in the next subsection. 



we can obtain the value of </> en d- We use this in the ex- 
pression for the number of e-foldings in the high-energy 
limit @ 



N 



16tt 2 

'3M B 8 ,„, 



„d y2 



(48) 



Taking N = 50, this gives the value of 4>5o and substitut- 
ing in the equations for ev and rjy one gets 



C. Specific examples 

To determine the typical size of the next-order correc- 
tion, we study the monomial potentials V oc <j) a for a = 2, 
4 and 6, assuming inflation takes place well within the 
high-energy regime. For comparison, we also calculate 
the magnitude of the correction for the standard cosmol- 
ogy 

To calculate the size of the correction term, we can 
use the slow- roll approximation for eH and 7/h, since any 
corrections to that will be of higher-order. The simplest 
approach is to rewrite ch and r]n in terms of the potential 
and its derivatives. We will make use of the equations 



H 



47T 

3M? 



V: 



5 

V; 
Ml H' 



(41) 
(42) 
(43) 



The first and second of these use the slow-roll approxi- 
mation, and the first and third use the high-energy ap- 
proximation. These enable us to obtain the relations 



ev,50 

W,50 



100 + 51a 

a — 1 
100 + 51a 



(49) 
(50) 



In the case of standard cosmology, the same calculation 
is carried out using the corresponding expressions for the 
slow-roll parameters and the power spectrum given by 
Refs. HQ 



V]i 2 {k) = [l-(2C+l)e H + C77H] 



2H 2 



Mi\H'\ 



(51) 



k=aH 



with en and rjn, defined as in Eqs. © and QJ, being the 
Hubble slow-roll parameters in the standard cosmology. 
To rewrite them in terms of the potential, we use the 
equations 0] 



H 2 
3H<p 



8tt 
-V: 



V; 



M 2 . 
— -H' . 

47T 



(52) 
(53) 
(54) 



e H ^ e v 



where the potential slow-roll parameters are 



(44) 

which arc the standard cosmology equivalents of 

same relations as in 

the previous case 



ev 



3M| V' 2 
16^2 ~y3 



3M| V" 
16^2 y 2 ' 



(45) 



being the high-energy limit of the parameters as defined 
by Maartens et al. |(| ■ Using them, we can write Eq. 1401) 
as 



7^ /2 = [l_(3C+l)e v + C77v] 



2H 3 



M§\H'\ 



(46) 



k=aH 



where we aim to calculate the square-bracketted term. 

We assume that observable scales crossed the Hubble 
radius 50 e-foldings before the end of inflation, and need 
to compute the slow-roll parameters at that time. We 
take the potential as 



V{(f>) = m(f> a , 



(47) 



where m is a constant. Putting this in Eq. (|45|l and set- 
ting ey = 1, which corresponds to the end of inflation, 



e H ^ e v 



with 



ev 



Mi V' 2 



'/v 



M 2 V" 



16tt V 2 ' 8tt V 

Then the power spectrum can be written as 

V][\k) = [l - (3C + l)ey + C m ] 



(55) 



(56) 



(57) 



k=aH 



where now N ~ -(8tt/M 4 2 ) rj™ d (V/V')d<f>. 

The results for both regimes are shown in Table [i] 
They show that the magnitude of the correction is similar 
in both cases, though it differs in detail. Nevertheless, 
these results confirm that the amplitude of the power 
spectrum is not changed significantly with respect to the 
slow-roll result by the higher-order correction. 
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TABLE I: Fractional correction to the power spectrum. 



a = 2 
a = 4 
a — 6 



standard cosmology 
ev.so Tv,50 correction 



ToT 1ST i- 0046 



51 



3 
102 



1.0019 



_2_ _2_ fl QQQ9 
103 103 u,roM 



high-energy braneworld 
ev,50 ??v,50 correction 



i 

101 



76 



3 
_>n:i 



1 

_MI2 



3 
31 1 J 



1.0082 
1.0085 
1.0086 



D. Correction from the high-energy approximation 

We end by mentioning that, especially when it is small, 
the higher-order slow-roll correction may well be sub- 
dominant to corrections coming from the high-energy ap- 
proximation not being exact. The slow-roll approxima- 
tion allows us to determine the size of such corrections; 
using Eq. I|17f) which is valid in any regime, we find the 
slow-roll perturbation amplitude to be 



,1/2 

n 



i + l (A 

2 \V 



2H 3 



Ml\H> 



(58) 



k=aH 



where we have expanded the result in terms of the pa- 
rameter X/V which is small in the high-energy regime. 
Because the term linear in X/V happens to cancel, the 
high-energy approximation is a better one than might 
have been expected, and we only need V > 10A to bring 
the correction from that approximation within one per- 
cent. More generally, comparison of Eqs. H4Q(I and l|58|l 
allows a test of when the higher-order slow-roll correction 
dominates the correction to the high-energy approxima- 
tion. 

If one were very ambitious, one could also attempt 
to generalize our higher-order slow-roll calculation to be 
valid in any regime, but given the limited ability of obser- 
vations to probe or distinguish amongst such small cor- 
rections, such a calculation does not seem worthwhile. 
The same remark applies to an attempt to calculate the 
higher-order gravitational wave spectrum correction; as 
gravitational waves are known to be subdominant such 
corrections are even less relevant, and also much harder 
to calculate due to the gravitational waves' ability to pen- 
etrate the bulk dimension. 



IV. CONCLUSIONS 



We have devised a Hubble slow-roll formalism for infla- 
tion in the RS-II braneworld cosmology, extending work 
by Hawkins and Lidsey 0, to define parameters 
which share the nice properties of those used in the stan- 
dard cosmology, which are recovered in that limit. As 
an application, we have computed the density pertur- 
bation spectrum in the high-energy limit, both exactly 
for power-law inflation and to higher-order for general 
slow-roll inflation models. To do so we have used the 
Mukhanov equation; while no one has yet been able to 
prove that this equation is still valid in the braneworld 
context, we have provided some evidence supporting its 
use. We have also quantified how well the high-energy 
approximation must hold in order for the higher-order 
slow-roll correction to be the dominant one. 

It is interesting to note that, having defined the slow- 
roll parameters en and ryn so as to give the usual spec- 
tral index formula for slow-roll perturbations, it turns 
out that the next-order correction is of the same form 
as in the standard cosmology. We are not aware of a 
physical reason which leads to this result. Neverthe- 
less, for a given choice of potential one expects that 
observable perturbations are generated at a different lo- 
cation on that potential depending on the braneworld 
regime, and so predictions for both the spectral indices 
(see e.g. Ref. [2(j) and for the higher-order corrections 
will be different. We have examined the magnitude of the 
correction for some simple potentials, and we conclude 
that there is no reason to believe that the higher-order 
correction might be more important in the high-energy 
regime than in the standard cosmology. As recent obser- 
vations including WMAP have restricted viable inflation 
models to regions close to the slow- roll limit, such cor- 
rections are expected to be small. 
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